Quantum dynamics of impurities in a ID Bose gas 
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Using a species-selective dipole potential, we create initially localized impurities and investigate 
their interactions with a majority species of bosonic atoms in a one-dimensional configuration during 
expansion. We find an interaction-dependent amplitude reduction of the oscillation of the impurities' 
size with no measurable frequency shift, and study it as a function of the interaction strength. We 
discuss possible theoretical interpretations of the data. We compare, in particular, with a polaronic 
mass shift model derived following Feynman variational approach. 
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I. INTRODUCTION 

Low-dimensional and strongly interacting systems 
have sparked intense scientific interest in recent years 
in diverse research fields, such as solid-state physics, 
nanoscience and atomic physics. In particular, strongly 
correlated systems display quantum phases dominated by 
quantum fiuctuations, such as the Mott-insulator phase 
[l[ and magncticlike ordered phases 0. In low dimen- 
sions, the interplay between interactions among particles 
and confining potential can enhance the effect of quantum 
correlations, resulting in peculiar regimes, such as Tonks- 
Girardeau or the sine-Gordon Q regime. Due to the un- 
precedented control over the interatomic interactions, the 
external trapping potentials, and the internal states of 
the atoms, ultracold atomic systems represent a versatile 
tool to explore these novel phenomena, and have already 
rovided the way to realize some of these quantum phases 
. A particular strength of cold atoms is the realization 
of systems that are hard to obtain in condensed-matter 
physics, such as multicomponent bosons of either differ- 
ent hyperfine states or different species 0, as well as 
the study of real-time dynamics of quantum many-body 
systems. 

In particular, cold atoms provide a way to realize a re- 
cently proposed new universal class of quantum systems, 
the so-called ferromagnetic liquids Q, occurring as the 
ground state of repulsively interacting two-component 
bosons. The time evolution of a spin flip, i.e., the ex- 
citation obtained by flipping one single bosonic "spin" 
from the fully spin-polarized ground state, thereby cre- 
ating a localized "impurity", leads to a class of dynamics 
completely different to that of the standard (Luttinger- 



liquid) theory but possessing remarkable universal fea- 
tures |8l4lC|. The physics of such impurities propagating 
in a sea of the majority species is further connected to 
several longstanding questions in strongly correlated sys- 
tems, both fermionic and bosonic. For fermionic systems 
and an immobile impurity, one example is the well-known 
x-ray edge problem, with the impurity leading to the An- 
derson orthogonality catastrophe[ll|. The motion of the 



impurity then strongly affects this physics, as was probed 
both for impurities in ^He [l3 and, more recently, for mi- 
nority species in interacting fermionic systems [isl , in 
the context of Fermi-liquid theory. Inextricably linked to 
the study of the impurity motion is the notion of po- 
larons [l^, i.e., the occurrence of density fluctuations of 
the majority species, both for fermions and bosons, lead- 
ing to renormalization of the impurity parameters, such 
as its mass. In one dimension, however, the nature of 
the bath is quite special, leading to the above-mentioned 
novel physical effects and, in particular, to subdiffusion at 
zero temperature. This problem is also directly related 
to the motion of a driven impurity which was investi- 
gated recently both theoretically |la - [l8| and experimen- 
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FIG. 1. (Color online) Ultracold Rb atoms with K impurities 
are loaded into an array of ID systems, "tubes" (left). The 
SSDP light blade spatially localizes the impurities into the 
center of the Rb tubes (right). 
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tally [19| . A system where the impurity-bath interactions 
are adjustable over a wide range paves the way for the 
study of a wealth of physical phenomena. 

In this work, we realize such a swtem using a species- 
selective dipole potential (SSDP) [20]. A minority species 
(K atoms) can diffuse into a majority species of (Rb) 
bosonic atoms, in a one-dimensional (ID) configuration 
(see Fig. [T]). We study how the interspecies interactions 
reduce the oscillation amplitude of the impurities' size 
"■(i) = \J (x^)- We compare the data in the light of a 
model based on polaronic mass shifts. 

The remainder of the paper is organized as follows: in 
Sec. II we describe the experimental procedure to pre- 
pare the one-dimensional (ID) systems of K impurities 
initially localized in the surrounding Rb bath and mea- 
sure their subsequent expansion; in Sec. Ill we report 
the experimental data; Section IV illustrates our theo- 
retical model based on a quantum Langevin equation. 
By means of the Feynman variational approach |2l| , we 
calculate the impurities' mass renormalization due to the 
interactions with the Rb bath. In Sec. V we compare the 
data with the results of the theoretical model, discussing 
the limitations and approximations involved. Finally, we 
summarize in Sec. VI. 



II. EXPERIMENTAL PROCEDURE 

The mixture of ^^Rb and ^-^K is first cooled to 1.5 /iK 
by microwave evaporation of Rb and sympathetic cooling 
of K in a magnetic trap, then loaded into a crossed dipole 
trap created by two orthogonal laser beams (A=1064nm, 
waists ~ 70 /im). Both species are spin polarized in 
their toj;^=:1) hyperfine states, featuring magnet- 

ically tunable interspecies interactions in the — 10 mT 
range. The mixture is further cooled by optical evap- 
oration performed in a uniform magnetic field of 7.73 
mT to adjust the interspecies scattering length a to the 
convenient value of 240 Bohr radii, ensuring both fast 
thermalization and a low rate of inelastic collisions. 

At this point we set the magnetic field to 7.15 mT, 
corresponding to vanishing interspecies interactions, and 
we adiabatically raise a vertical (z) ID optical lattice 
(waist 170 /xm) to 15 (6.5) E^cc for Rb (K) in 200 ms. 
The lattice wavelength A=1064 nm results in a 2.3 times 
smaller lattice height for ^^K than for ®^Rb in units of 
the respective recoil energies (ii'rec = h? /2m\^). After 
adiabatic extinction of the dipole trap, both species move 
downward as collisions disrupt the gravity-induced Bloch 
oscillations [l^ but, due to largely different tunneling 
times, K drops faster and, by adjusting the fall time, we 
reduce the gravitational sag between the two species in 
the dipole trap. 

Subsequently, we further raise the vertical lattice to 
60 E'roc for Rb and adiabatically switch on an additional 
standing wave with equal waist and strength along the y 
direction. Thus we create the array of ID systems (tubes) 
by means of a two-dimensional (2D) lattice. The lat- 
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FIG. 2. (Color online) Time sequence of the experimental 
procedure used to prepare the K impurities in the middle of 
ID tubes filled with Rb atoms. 

tice transverse harmonic oscillator frequency, u)^/2t: =34 
(45) kHz for Rb(K), exceeds both the temperature and 
the chemical potential, thus ensuring the ID regime for 
both species. 

The residual trapping frequency along the direction x 
of the tubes, 62 (87) Hz for Rb(K), is due to the inhomo- 
geneous transverse profile of the lattice beams. An SSDP 
elliptic beam, orthogonal to the tubes, is then turned on 
in 50 ms, compressing K at the center of Rb and leav- 
ing the latter nearly unperturbed. This light blade has 
waists of 15 and 75 /im {x and z direction, respectively), 
wavelength of 770.4 nm, and power equal to 0.6 mW. 
Correspondingly, the depth and frequency are 11 /xK and 
a;^/(27r) = 1.0 kHz. 

Finally, the interaction strength 5iD,KRb between K 
and Rb, is brought to the desired value by linearly 
ramping the magnetic field to its final value Bd- In- 
stead, the interaction strength of the Rb atoms .giD.Rb = 
2.36 X 10^"^^ Jm is independent of the applied magnetic 
field Bd. 

The impurities' dynamics is initiated by rapidly extin- 
guishing the light blade with a linear ramp of 0.5 ms. 
Images of the impurities are taken once their motion has 
been frozen by suddenly adding a tight optical lattice 
along the tubes axis allowing for magnetic field to be ex- 
tinguished in 15 ms and for atoms to be repumped in the 
hyperfine level suitable for imaging without significantly 
affecting their density distribution. The full experimental 
procedure is shown in Fig. [5] 

Typically, we measure (1.8 ±0.2) x 10^ ^'^Rb and (6 ± 
2) X 10^ ^^K atoms at 140 nK before loading the tubes. 
From the analysis of time-of-flight images, we obtain a 
temperature of the Rb sample in the ID tubes of 350(50) 
nK. We estimate the Rb peak filling and density to be 180 
atoms/tube and niD,Rb — 7/im~^, with a filling-averaged 
Lieb-Liniger parameter (mijf,5iD,Rb)/(^^'T-iD,Rb) — 1; the 
peak K filling is approximately 1.4 atom/tube. 



III. EXPERIMENTAL DATA 

We record the oscillation of the impurities axial size 
a{t) ~ (a;2) along the tubes as a function of time, 
through resonant in situ absorption imaging. We fit the 
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FIG. 3. (Color online) Upper panel; oscillations of the K impurities axial width a — \J {x^) after extinction of the confining 
SSDP, for different interaction strengths with the surrounding Rb bulk, 77 = .giD,KRb/<7iD,Rb- Solid lines represent a fit to 
data, triangles the Rb axial size. Inset images show the observed in-situ density distributions of Rb (red, left) and K (blue, 
right) within a window of 200^m xl50/im. Lower panel: values of the fit parameters uj/2ti, 7 and /3 (left to right, see text for 
definition). 



absorption profiles by a two-dimensional Gaussian func- 
tion. Figure [3] shows the evolution of the impurities axial 
size a it) for four different values of the ID interspecies in- 
teraction strength gio.KRb- Given the magnetic field -B^, 
hence the interspecies scattering length a [23|, gio.KRb is 
calculated following [U and expressed in units of jfiD.Rb 
as .9iD,KRb = '75iD,Rb- We find the oscillation amplitude 
of (T{t) to be large for 77 ~ (no interactions between 
impurities and ID gas) and to decrease with 77. Remark- 
ably, for the largest oscillation strength, the impurities 
oscillations are always confined inside the bath. In ad- 
dition, we observe that, in addition to oscillating, a{t) 
increases linearly (at least initially) over time. In order 
to extract quantitative information, we fit the experimen- 
tal data with a damped sine function with linear baseline: 
cr(t)=(7i + pt - Ae-^"^*- cos(v/l - I'^uj {t - to))- Since in 
the following we will mainly focus on the oscillation am- 
plitude, in Fig.|3]wc show the fit parameters, aj/27r, 7, /?. 

We first notice that the oscillation frequency does not 
shift, within our error bars, with gio.KRb: such a surpris- 
ing feature is further discussed below. Then we observe 
that, even in the absence of interactions at rj = 0, a 
residual damping occurs, likely due to intraspecies col- 
lisions in tubes with more than one K atom. We also 
notice that, during the oscillations, the impurity cloud 
can grow larger than the Rb sample, thus exploring the 
Rb inhomogeneous density profile. Such inhomogeneity 
along the tubes axis adds to the inhomogeneity of the 
tube's filling. In the following, we account for the inho- 
mogeneity using the local density approximation, i.e., by 
considering homogeneous systems with different values of 



Rb linear density. As K impurities are initially confined 
in the middle of Rb, the first oscillation is less affected by 
Rb inhomogeneous density than the following. Therefore 
we focus our analysis on the value of the maximum size 
reached in the first oscillation, after 3 ms of expansion: 
CTp = (T{t = 3 ms). Such value closely refiects the oscilla- 
tion amplitude A as, for short times, both the slope and 
the damping have negligible effect. 
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FIG. 4. (Color online) Experiment: impurities' axial size 
at the first oscillation maximum ap, normalized so that 
ap — 1 for rj = 0, versus the coupling strength parameter 
77 = giD.KRb/ffiD.Rb, for attractive (left, 77 < 0) and repul- 
sive (right, »7 > 0) interactions (circles). Theory: rriK /rn*j^ 
calculated with impurity-bath coupling giD.KRb derived from 
two-body scattering [2J] (solid line) and with gio.KRb multi- 
plied by a fit parameter (dashed line). Both curves are com- 
puted for a Rb density of 7 atoms/^m and 7 = 0. 
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FIG. 5. (Color online) Amplitude of the first oscillation of the 
axial size of K impurities (right) following different prepara- 
tion sequences (left): (a) 77 always negligible; (b) 77 large while 
K in the SSDP light blade and negligible during expansion at 
t > 0; and (c) 77 always large. 



Figure S] shows the dependence of ap on the relative 
coupling strength rj: values of 77 of the order of the unity 
are sufBcient to induce a significant amplitude reduction, 
if compared to the non interacting case 77 = 0. It is impor- 
tant to remark that, if mean- field interactions dominated, 
we would observe an opposite behavior between positive 
and negative rj values. As ap decreases with |r7|, indepen- 
dently of its sign, we conclude that mere mean-field in- 
teractions cannot explain our experimental findings. We 
also notice a saturation behavior for 77 > 4 that cannot be 
due to the washing out of the Feshbach resonance caused 
by magnetic field instability (~ 10 /xT). Instead, it might 
be related to a crossover to the three-dimensional (3D) 
regime, as we estimate that for 77 = 15 the mean-field 
interaction energy of the K impurities, (7iD,KRb'^iD.Rb, 
equals the band-gap of the 2D lattice. As our theory 
cannot describe this ID to 3D crossover, and furthermore 
the inclusion of the ID bound states between K and Rb 
for 77 < al so cannot be adequately captured, we have 
calculated rriK / only in the range < 77 < 8 in Fig. 
3. 

To some extent, the amplitude of the first oscillation 
could reflect different preparation temperatures of the K 
sample in the light blade, as the efficiency of the thermal- 
ization with the Rb background depends on the strength 
of interactions. Thus, we recorded the first oscillation 
of K impurities prepared at large 77 and expanding at 
zero interactions, and compared with the oscillations of 
samples prepared and expanded at large (or zero) inter- 
actions. We notice from Fig. [S] that the preparation has 
a nearly negligible impact on the oscillation amplitude, 
insufficient to explain the large differences observed be- 
tween high and small 77. Moreover, incomplete thermal- 
ization cannot explain why at large interactions Up lies 
below the Rb axial width, as shown in Fig. [31 

To summarize, the key experimental findings are (a) 
even for large 77 the oscillation frequency of a{t) does not 
deviate measurably from that of the noninteracting case 
(which is 2a;); (b) the larger I77I, the lower the amplitude 
of oscillation initially; (c) in the interacting case, besides 
oscillating, the impurities width increases over time in 
a seemingly linear fashion; and (d) at long times a{t) 



equilibrates to about the same value, independently of 77. 



IV. THEORETICAL MODEL 

A complete explanation of the observed phenomena is 
an interesting and open problem. Motivated by the ob- 
served shape of cr(i), we give here a semiempirical anal- 
ysis of the oscillation through the model of a damped, 
quantum harmonic oscillator in contact with a thermal 
bath, i.e., a quantum Langcvin equation (25j . x{t) = 
p{t)/m*j^, p{t) = —k*x{t) — jp{t) + £,{t), as a framework 
to order and provide a first interpretation of the find- 
ings. Here, x{t) and p{t) are respectively the position and 
momentum operators of one ^^K atom whose mass nij^ 
and spring constant k* are renormalized by interactions 
with the bath. In addition to the exponential damping 
—jp(t), we account for the fluctuations of the bath by 
the noise operator (,{t), whose correlator {£,{t)£,{t')) cx 7 
is approximated by using the thermal quantum statistics 
of a set of harmonic oscillators at constant spectral den- 
sity [25|. The operator £^{t) causes the gradual increase 
of the impurities' width, corresponding to heating from 
the thermal component of the bath. 

Various effects could account for the change of am- 
plitude. A dominant mean-field interaction for 77 > 
would make the effective potential for the K shallower, 
resulting in ap increasing with rj, at odds with observed 
behavior. Another possibility is the damping from the 7 
term. However this is also inconsistent with the observed 
data: within the quantum Langevin approximation, and 
at the experimental temperatures, any but very small val- 
ues of 7 result in strong heating contributions from the 
correlator {£,{t)(,{t')), such that overall ap increases with 
77. 

At variance with the two above-mentioned effects, a 
mass renormalization of the K impurities due to polaronic 
effects [21| is a very good starting point to explain the 



decrease of (T„. 



Polaronic mass shift 



We calculate the mass shift M 



Tij^ — tuk usmg 

Feynman's variational theory for the polaron [2l| . Taking 
the combined impurity-bath Hamiltonian 



H 



2mK 



ikx 



bk + b 



(1) 



where the Rb bath has been approximated as a 
Tomonaga-Luttinger liquid with linearized density fluc- 
tuations around a homogeneous background density, de- 
scribed by the bosonic operators bk bl with = 

Vk = .giD,KRb(-ft:|A:|/27rL)^/^e-l'=l/2'=<=, and where Vs is 
the sound velocity of Rb, K the Luttinger parameter, re- 
lated to the effective Rb-Rb interaction and kc the cutoff 
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FIG. 6. (Color online) Square root of the inverse impurities' 
renormalized mass versus interspecies interaction strength 77, 
calculated for different values of Rb density. This quantity 
directly compares with the measured K width (jp. 



momentum. These quantities can be related to a homo- 
geneous ID Rb gas with arbitrary contact interactions 
[26|; as a working theory for inhomogeneous ID systems 
is difficult to obtain, we use this homogeneous theory in- 
stead, evaluating it for different Rb densities in the center 
of the trap. We find that the results are not that strongly 
depending on the precise value of the Rb density. 

Evaluating the full partition function Z = tr (^e^^^ 

as a path integral and integrating out the Rb bath yields 
an action S for the K atom that is non-local in time: 



S = l dr^i\r)- 



(2) 



n\ ik[x(T)-x{T'}] 



K = exp 



2{mK + M) 
cosh{Wahf3/2) 



M 

np rriK 
cosh[Wa{ni3/2 - u)] 



Wa sm\i{hpWa/2) 



(4) 



which we then minimize numerically to obtain the mass 
shift M. 

We apply this theory to an experiment with an inho- 
mogeneous density distribution of bath particles (due to 
the external confinement of the Rb atoms) . For this ap- 
proach to be reasonable, the mass shift should not depend 
too strongly on a shift in Rb density, which Fig. |6] shows 
to be indeed the case (all other parameters are as in the 
experiment). 

Note that one can easily add a term corresponding to a 
parabolic confining potential in the action ([2]) and repeat 
the variational procedure to find the effective mass. Quite 
intuitively, one finds that the tighter the confinement, the 
smaller is the mass renormalization. This strongly sug- 
gests that before the blade is released one should consider 
that the K impurity has its bare mass. This is clearly a 
point that would need to be more quantitatively tested 
by, e.g., a Monte Carlo calculation. 



COMPARISON OF DATA AND THEORY 



where G{u, k) = cosh(efc(M-/i/3/2))/ sinh(;i/3efc/2) The 
mass shift enters through the Feynman trial action 6*0, 



^0^ dT——x (7 



8 



sinh(VF/3ft/2) 



which depends on two parameters, the mass shift M and 
W . which are chosen by minimizing the trial free energy 
Fo + {S - So) /{hp) >F = -(l//?)lnZ. 

The explicit expression of the trial free energy is 



Strongly simplifying the bath (see [25|) for a small 7 as 
the data suggest we obtain the functional form for <T{t) 
used to fit the oscillations in Fig. [31 The theoretically pre- 



dicted decrease of ap is proportional to 1/ , 



There- 



fore, in Fig. 2] we compare the measured <7p normalized to 
the value at zero interactions to the calculated values of 
\J mK Ittl*^: the theoretical curve reproduces the experi- 
mental trend and magnitude of the amplitude reduction 
fairly well. 

As the variational approach cannot take the inhomo- 
geneity of the bath into account, we approximate the 
translationally invariant theory with representative val- 
ues for the Rb density in the center of the trap. We find 
that, in the parameter range of the experiment, the mass 
shift is rather insensitive to the Rb density. 
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We also considered the possibility that (7iD,KRb in the 
experiment {i.e., at finite density) may be different from 
the one obtained in the standard ID two-body scatter- 
ing theory [23, H^] ■ While those formulas have indeed 
been demonstrated to well describe the position of the 
confined-induced resonances (CIR) also in the regime of 
finite density, it is currently unknown how accurate it is 
for the actual value of gio.KRb away from the CIR. We 
find that a least square fit of the theory to data with a 
value of 3.15(7iD.KRb fits the experiment better. 

Concerning the mass shift theory for cTp, we must note 
that for the trap parameters and temperatures it relies 
on the assumption that the mass of K is not or is only 
insignificantly renormalized inside the light blade dur- 
ing preparation. This is supported by calculations of the 
mass rcnormalization that include a parabolic trap for 
the K, which show the mass shift becoming smaller the 
deeper the trap, eventually becoming zero. This simple 
extension of Feynman's approach is, however, ill suited to 
quantitative calculations of concrete values of the mass 
shift in a tight confinement, and we therefore must leave 
the quantitative test of this assumption to strong numer- 
ical methods. 

The above interpretation based on the polaronic mass 
shift also predicts an increase of the oscillation period, 
which for large 77 would be beyond the 5% experimen- 
tal error bars. To reconcile this apparent contradiction 
with the data, we need to consider the physical processes 
resulting in an effective upward rcnormalization of the 
spring constant k* . Such processes are more easily visu- 



alized in the limit of impenetrability: when the impuri- 
ties cannot get across the bath atoms, their displacement 
from the trap bottom costs the additional harmonic po- 
tential energy of the bath atoms moving uphill. While 
this is observed in proof-of-principle simulations for sim- 
plified models, quantitatively estimating the value of k* 
is difficult with current methods and will require further, 
more refined analysis. 

VI. CONCLUSIONS 

In conclusion, we have investigated the dynamics of 
impurities in ID atomic samples as we varied their in- 
teraction strength with the surrounding bath of weakly 
interacting bosons by more than two orders of magnitude. 
We observe the amplitude of the impurities' quadrupole 
oscillations decreasing with the absolute value of the in- 
teraction strength, a fact unexplained by mere mean-field 
interactions. An analysis of the data in the light of a 
simplified quantum Langevin equation suggests that po- 
laronic mass shift plays a role in the amplitude reduction. 
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